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A compact steep spectrum radio source in NGC 1977 
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ABSTRACT 

Context. A compact steep spectrum radio source (J0535-0452) is located in the sky coincident with a bright optical rim in the Hn 
region NGC 1977. J0535-0452 is observed to be < 100 mas in angular size at 8.44 GHz. The spectrum for the radio source is steep 
and straight with a spectral index of -1.3 between 330 and 8440 MHz. No 2 /im IR counter part for the source is detected. These 
characteristics indicate that the source may be either a rare high redshift radio galaxy or a millisecond pulsar (MSP). 
Aims. To investigate whether the steep spectrum source is a millisecond pulsar. The optical rim is believed to be the interface between 
the Hn region and the adjacent molecular cloud. If the compact source is a millisecond pulsar, it would have eluded detection in 
previous pulsar surveys because of the extreme scattering due to the H n region-molecular cloud interface. 

Methods. The limits obtained on the angular broadening along with the distance to the scattering screen are used to estimate the pulse 
broadening. The pulse broadening is shown to be less than a few msec at frequencies J 5 GHz. We therefore searched for pulsed 
emission from J0535-0452 at 14.8 and 4.8 GHz with the Green Bank Telescope (GBT). 

Results. No pulsed emission is detected to 55 and 30 jiiy level at 4.8 and 14.8 GHz. Based on the parameter space explored by our 
pulsar search algorithm, we conclude that, if J0535-0452 is a pulsar, then it could only be a binary MSP of orbital period 5 5 hrs. 
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1. Introduction 

NGC 1977 is an Hn region located at the northern edge of the 
Orion molecular cloud at a distance of about 0.5 kpc from the 
Sun. The region has an emission measure of 5 x 10 3 pc cm -6 
and is believed to be an example of an interface between the 
Hn region and the adjacent molecular cloud (Shaver & Goss 
1970). A 6-cm VLA image of the H n region made by Kutner et 
al. (1985) shows a strong peak that coincides with the brightest 
part of the optical bright rim. They investigated whether the ra- 
dio peak is due to a compact H n region or dense clump of gas 
ionized by the central star in NGC 1977. Both these possibilities 
were ruled out based on the expected radio and IR emission from 
such ionized gas. Subrahmanyan, Goss, & Malin (2001) imaged 
NGC 1977 with the VLA at 330, 1420, 4860 and 8440 MHz. 
They detected the bright radio peak as an unresolved object des- 
ignated as J0535-0452. The estimated flux density of the com- 
pact object is 320, 50, 1 1.4 and 5.2 mJy at 330, 1420, 4860 and 
8440 MHz respectively. The spectrum of the object is straight, 
with no turn over down to 330 MHz. The spectral index is -1.3 
over this range of frequencies. The 8440 MHz observations are 
used to obtain an upper-limit on the source size of ~ 100 mas. 
The implied brightness temperature is 7 x 10 s K. Subrahmanyan 
et al. (2001) have concluded that the object is a non-thermal 
source. 

Kutner et al. (1985) argued that the non-thermal source 
could be an extragalactic background source located behind 
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NGC 1977. Compact steep spectrum extragalactic sources with 
spectral index close to -1.3 are typically high redshift objects 
(> 2; Miley & De Breuck 2008). For these redshifts, the size of 
J0535-0452 is J 780 pc, assuming standard cosmological pa- 
rameters. Compact steep spectrum sources of this size typically 
show spectral turn over near a few GHz, which is not the case 
for J0535-0452. Thus if J0535-0452 is an extragalactic source, 
then it should be one of the rare steep-spectrum, high redshift 
radio galaxies (for example, see Kaplan et al. 2000). 

Another possibility is that the compact source in NGC 1977 
is a pulsar. The spectral indices of both normal and millisecond 
pulsar are in the range -1 to -2.5 with a mean value of -1.65 
(Karmer et al. 1998). Thus the compactness and spectral index 
of J0535-0452 may suggest that it is a pulsar. An earlier attempt 
to detect pulsed emission near 330 MHz was not successful 
(Ramachandran, R. 1995; unpublished). However, J0535-0452 
would have eluded detection in earlier pulsar searches due to 
extreme scattering. We show in Section [2] that the best frequen- 
cies to search for pulsed emission from this pulsar are > 5 GHz 
depending on the pulsar period. We used the GBT to search 
for pulsed emission at frequencies near 4.8 and 14.8 GHz, but 
did not detect pulsed emission. The observations and results are 
given in Section[3]and[4] respectively. 



2. On observing pulsed emission from J0535-0452 

An attempt was made in 2003 to image J0535-0452 using the 
VLBA at 8.4 GHz (Roshi et al. 2003 unpublished). The source 
was not detected in this observation. A possible interpretation of 
the non-detection in the VLBA observations is that the source 
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Fig. 1. (a)Schematic showing the scattering geometry. The H n region-molecular cloud interface (scattering screen) is located at 
~ 450 pc from the observer and the distance to the compact source is > 2.5 kpc. (b) The expected pulse broadening for three 
frequencies (330 MHz, 5 and 14 GHz) is shown in the pulsar period-angular broadening parameter space. The curves are the locus 
of points for pulse broadening equal to the pulsar period. For each frequency the parameter space above the curve can be probed by 
the pulsar search observation. The pulse broadening estimation includes contribution from an assumed intrinsic duty cycle of 10% 
of the pulsar period and broadening due to DM. 
For these calculation, the smallest pulsar period is taken as 1 msec. The range of angular broadening is constrained by observations between 10 
and 100 mas. For all frequencies 5 14 GHz, the pulse broadening is dominated by scattering. At 14 GHz, for angular broadening i 60 mas, the 

pulse broadening is limited by the assumed intrinsic duty cycle. 



would have angular broadened due to scattering in the H n re- 
gion - molecular cloud interface. Assuming that the source is 
scatter broadened, a lower limit can be obtained for angular 
broadening which turns out to be ~ 10 mas. If J0535-0452 is 
a pulsar, then the pulsed emission will be broadened due to scat- 
tering in the Hn region - molecular cloud interface. We esti- 
mate the pulse broadening using the measured limits on angu- 
lar broadening of J0535-0452 along with the source-scattering 
screen geometry shown in Fig QJ. The distance to the compact 
object, Dp U i, as measured from H i absorption studies is > 2.5 
kpc (Subrahmanyan et al. 2001) and the distance to NGC 1977, 
D m, is ~ 450 pc (Genzel et al. 1981). The e~ l pulse broaden- 
ing time in sec due to scattering is obtained using the equation 
(Cordes & Lazio 1997) 

t _ D hii D P«l liAGm JfcHzY 4 

s D pul -D htt 81n(2)c v 8.44j ' 

where c is the velocity of light in cm s , the units of distances 
are in cm, fcHz is the observing frequency in GHz and 8 s $achz is 
the angular broadening at 8.4 GHz in rad. The pulse broadening 
time is considered to scale with the 4'' 1 power of frequency, ap- 
propriate for strong scattering. 9 s ,?,aghz range between 10 and 
100 mas, which are obtained from the VLA (Subrahmanyan 
et al. 2001) and VLBA observations near 8.4 GHz (Roshi et 
al. 2003 unpublished). In Fig QJ) we plot the estimated pulse 
broadening time for three observing frequencies (330 MHz, 5 
and 14 GHz) in the pulse period- angular broadening parameter 
space. The pulse broadening time in this plot also includes an 
assumed intrinsic duty cycle of 10% of the pulsar period (mini- 
mum pulsar period is taken as 1 msec). The spectral resolution 
used for the calculation is 0.78 MHz. For a given dispersion 
measure (DM), this finite spectral resolution produces temporal 



broadening, which is also added to the plotted pulse broadening. 
Any temporal broadening due to DM griding in the pulsar search 
algorithm is neglected. 

Continuum observations toward J0535-0452 can be used to 
estimate limits on the DM. The measured spectrum of the com- 
pact source does not show any turnover at 330 MHz due to the 
continuum optical depth of the H n-region-molecular cloud in- 
terface. This fact is used to obtain an upper limit on the emis- 
sion measure (EM) of ~ 10 4 crrr 6 pc for an assumed ionized 
gas temperature of 8000 K and t c ~ 0.05 at 330 MHz. The 
3 sigma uncertainty of the 330 MHz map of NGC1977 region 
(Subrahmanyan et al. 2001) is used to infer an upper limit on r c 
of ~ 0.05. A DM corresponding to a given EM can be obtained 
as DM - y/EM (j)L, where <p is the filling factor of the ionized 
gas in the intervening H n region and L is the line of sight extent 
of the H n region. Taking <p ~ 1 and L ~ 10 pc, a good fraction 
of the size of the Orion molecular cloud, the DM we get is ~ 300 
pc cm 3 . Contribution to DM due to the distributed ionized gas 
has to be added to the above estimate. NE2001 model (Cordes & 
Lazio 2002) provides a total DM of ~ 50 pc ctrT 3 in the direction 
of NGC1977. Thus the upper limit obtained on DM is ~ 350 pc 
cm -3 . For estimating pulse broadening we used DM values up to 
650 pc crrr 3 . 

The result of the pulse broadening calculation shows that for 
all frequencies 5 14 GHz the pulse broadening is dominated by 
scattering (see Fig|TJ)). At 14 GHz, for 9 s ,%aghz < 60 mas, pulse 
broadening plotted in Fig QJ) is limited by the assumed intrinsic 
duty cycle. For each observing frequency indicated in Fig QJ) 
pulsed emission can be detected in the parameter space above 
the corresponding curve. 

The estimated continuum flux density at 14.0 GHz is 2.5 
mJy. If the observed continuum emission is the mean flux den- 
sity of the pulsar, then the pulsar can be detected at frequencies 
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5 and 14 GHz (depending on the pulsar period) with signal-to- 
noise ratio > 100 in 1 hr of observing time with the GBT. Thus, 
as seen in Fig.QJ), such an observation will be sensitive to both 
MSP and ordinary pulsars. Note that earlier searches would not 
have detected pulsars for the measured angular broadening limits 
(see Fig[Tp). 
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3. The GBT Observations and Data Analysis 

The GBT observations on J0535-0452 were made at 4.8 and 
14.8 GHz on 15 May 201 1. The GUPPI backend with 800 MHz 
bandwidth, 1.6 MHz spectral resolution and 41 /isec time res- 
olution was used for the observations. Flux density calibration 
was done using the source 3C161 for both frequencies. The flux 
densities of 3C161 were taken to be 6.7 and 2 Jy at 4.8 and 14.8 
GHz (Ott et al. 1994). The measured telescope gains were 2 and 
1.9 K/Jy and system temperatures were 18.5 and 26 K at 4.8 and 
14.8 GHz respectively. Reference pointing was done on B0540- 
0415 at 14.8 GHz and on 3C161 at 4.8 GHz. The total on-source 
observing time at 14.8 GHz is 27 minutes and that at 4.8 GHz is 
29 minutes. The data processing was performed using PRESTCQ 
(Ransom, Eikenberry & Middleditch 2002). A DM range of to 
720 pc crrT 3 was searched with full time resolution. The value 
of 720 pc cm~ 3 is about two times larger than the upper limit on 
DM estimated in Section [2] In addition we searched for pulsar 
acceleration in the range to 2. 1 xl0 4 x P m s -2 , where P is the 
period of the pulsar in sec. 



4. Results and Conclusion 

We did not detect pulsed emission from J0535-0452. The upper 
limit obtained for pulsed emission at 4.8 GHz is 55 fiJy for a 
pulsar period of 10 msec and pulse broadening of 7 msec cor- 
responding to an angular broadening of 60 mas (see FigQJ)). At 
14.8 GHz, the upper limit obtained is 30 //Jy for a pulsar period 
of 3 msec, typical for MSP, assuming 10% intrinsic duty cycle. 

Acceleration searches are most sensitive to binary pulsars 
with orbital period greater than 10 times the observation time 
(Ransom, Cordes & Eikenberry 2003). This means that our data 
analysis will rule out compact sources in binary system with or- 
bital period Z 5 hrs (assuming circular orbit). We have searched 
for accelerations up to 2.1 xl0 4 x P m s -2 . For orbital period of 
~ 5 hrs, this upper limit rules out companion object of mass 5 
3 M for an assumed pulsar period of 3 msec. If the companion 
objects is of higher mass (ie > 3 M ) then for orbital periods > 5 
hrs the object cannot be an ordinary star. This is because the or- 
bital radius becomes comparable to the size of the star. The com- 
panion object could be a stellar mass black hole. However, we 
are not aware of any other observational evidence for the pres- 
ence of a stellar mass black hole in the direction of NGC1977. 
Thus the only possibility that remains is that J0535-0452 is a 
binary MSP with short orbital period (ie J 5 hrs). 
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